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?N ■ Abstract 

We consider models with multiple Higgs scalar gauge singlets and the resulting restrictions on 

■ the parameters from precision electroweak measurements. In these models, the scalar singlets mix 
with the SU{2)i Higgs doublet, potentially leading to reduced couplings of the scalars to fermions 
and gauge bosons relative to the Standard Model Higgs boson couplings. Such models can make 
the Higgs sector difficult to explore at the LHC. We emphasize the new physics resulting from the 

■ addition of at least two scalar Higgs singlets. 
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I. INTRODUCTION 



One of the major goals of the Large Hadron Colhder is to probe the electroweak symmetry 
breaking sector. The simplest implementation of the symmetry breaking utilizes a single 
SU{2)l scalar Higgs doublet. In this minimal case, the couplings of the Higgs boson to 
fermions and gauge bosons are fixed in terms of the particle masses and the phenomenology 
has been extensively studied. It is of interest, however, to study extentions of the Higgs 
doublet model and to examine which possibilities are allowed by current data and how LHC 
Higgs phenomenology is affected. The most straightforward possibility for enlarging the 
Higgs sector is to add some arbitrary number of scalar singlets which couple only to the 
Higgs doublet. 

The phenomenology of models with one scalar singlet in addition to an SU{2)l doublet 
has been examined by many authors[ll, i, Si, H, Q. The case with one additional scalar 
is similar to that of models with a radion[7|. For a supersymmetric model, the addition 
of a gauge singlet scalar superfield leads to the NMSSMjsl, 0, [lo|, which solves the so- 
called "/i" problem of the MSSMfTT]. Alternatively, scalar singlets have been advocated 
as a signal for a hidden world which interacts only with the scalar sector of the Standard 

Model m M, Q El 



In this paper, we consider non-supersymmetric models with additional scalar gauge sin- 
glets. In the case where there is a Z2 symmetry in the scalar sector, this class of theory 
generically leads to a dark matter candidate, which is the lightest scalar singlet. Without a 
Z2 symmetry, the scalar singlets can mix with the Standard Model Higgs doublet and there 
is no dark matter candidate. It is this alternative which we consider here. The existence 
of multiple scalar singlets leads to changes in the scalar interactions with gauge bosons and 
fermions. Many authors have considered the case where the lightest scalar has a mass on 
the order of a few GeV and attempted to construct scenarios which evade the LEP direct 
Higgs production bounds^, 1^. We consider an alternative case where all the scalars are 



heavier than the LEP lower bound on the Standard Model Higgs, Mh,sm > 114 Ge\^|17l . 

The existence of scalars heavier than the LEP bound is restricted by electroweak precision 
measurements 18. [TqI . 20| . Since the dependence of the electroweak measurements on the 



scalar masses is logarithmic, it is possible to make quite significant changes in the scalar 
sector and still be consistent with precision dataHEi, [23j. We compare the predictions 



of models with multiple Higgs scalars with the restrictions obtained from the S, T, and U 



parameters [24|. We examine the cases with one and two scalar singlets and derive some 
general restrictions on the properties of models with extra scalar singlets. In the Standard 
Model, precision electroweak measurements restrict the Higgs mass to be less than about 
185 GeV, Mh,sm < 185 GeV^[l9[. We consider the possibility of discovering a Higgs-like 
boson with a mass significantly larger than allowed in the Standard Model in a theory with 
multiple scalars. As more and more singlets are added, the couplings of the individual 
scalars to the fermions and gauge bosons become weaker and weaker and heavy scalars can 
potentially be compatible with precision measurements. We are motivated by the analysis 
of Ref.[H which attempted to hide the Higgs signal at the LHC by introducing multiple 
Higgs scalars. This reference concluded that a model with three scalars with masses in the 
120 GeV region could evade discovery at the LHC^. 



^ Ref. [l| found that a model with three scalars with masses mo = 118 GeV, mi = 124 GeV, and 
m2 = 130 GeV would elude detection at the LHC with L = 100 fb^^. 
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In Section [ITl we summarize the class of models which we consider in this note, while 
Section IIIII contains our results for the S, T, and U parameters. Our results for one and 
two singlets are contained in Section llVt along with a discussion of the phenomenological 
implications of our results for Higgs searches at the LHC Technical details are summarized 
in two appendices. Section |V] contains some conclusions. 



II. THE MODELS 



We consider a class of models with N scalar singlets, Si, along with an SU{2)l doublet, 



H 



Si = Si + V, 



Si 



(1) 



We assume that the scalar potential is such that all scalars get a vacuum expectation value 
(VEV), 



(2) 



Since the singlets do not couple to the SU{2)l x f/(l)y gauge bosons, they do not contribute 
to My/ and Mz and hence vh must take the Standard Model value, vh = 246 GeV . The 
VEVs are determined from the scalar potential. 



|2\2 



I Si + —Ylijfi^jSiS. 



H 



Si + YjijAd^jSiSj + Yjiji^XijifSiSjSk + Yjiji^iXiji^iSiSjSkSi 



(3) 



Note that we make no assumptions about possible Z2 symmetries in the scalar sector and 
in general h will mix with the Sj scalars to form the mass eigenstates. 

The + 1 scalar mass eigenstates are defined to be = 0...N, with masses, rrii. 
We assume that mo is the lightest scalar. The mass eigenstates are related to the gauge 
eigenstates by an (A^ + 1) x [N + 1) unitary matrix V, 

Sl 



V 



(4) 



\4>N J J 

Our results are expressed in terms of the elements of the mixing matrix V, which can 
be calculated in any given model. The couplings of the scalars to the gauge bosons and 
fermions are^ 



■^^M ^ AT ^ 



H i,j=0,N 



(5) 



Note that the Goldstone bosons have Standard Model couphngs. 
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The production rates of the 0i are suppressed by | Voi p relative to the Standard Model 
Higgs boson production rates. The branching ratios of the lightest scalar, 0o, to Standard 
Model particles are identical to the Standard Model branching ratios, while the branching 
ratios for the heavier scalars depend on whether the channels (pi — > (pjcpk are kinematically 
accesiblej2|]. 



III. LIMITS FROM PRECISION ELECTROWEAK MEASUREMENTS 



The limits on the parameters of the scalar sector from precision electroweak mea- 
surements can be studied assuming that the dominant contributions resulting from the 
expanded scalar sector are to the gauge boson 2-point functions[2J, |25|, U'^yip'^) = 
IIxy{p'^)9'"' + Bxy{.P^)p^p'', with XY = 77, 7Z, ZZ and W+W- . We define the S,T and U 
functions following the notation of Peskin and Takeuchi 2J] , 



aS 



V Ml 



r-2 - e2 



cese 



UzziMl) -UzziO) -U^,iMl] 



n^z(M|)-n^^(o) 



aT 
aU 



UwwiO) n^z(O) 2se U^zjO) 
■ ce Ml 



Ml 



Iiww{Ml,)-Yiww{^) 



UzziMl) - UzziO) 



Ml 



^ , U,z{Ml) - U,z{0) ^ „2n,,(M|) 



Ml 



Ml 



(6) 



where sg = sin 9^ and cg = cos 9\y and any definition of sg can be used in Eq. [6] since the 
scheme dependence is higher order. 

The scalar contributions to S, T, and U from loops containing the (pi are gauge 
invariant [2^ and can be found from Appendix 1 of Ref. [27|] or from Ref. [28]^. 

= -S, I Voi P (So(0,m„Mz) -5o(Mz,m,,Mz) 

B22{Mz,mi,Mz) ~ B22{0,mi,Mz) 



Ml 



^ I Voi P |-5o(0, mi, Mw) + \bo{0, m„ Mz) 



1 



^22(0, rrii, Mw) - ^22(0, rrii, Mz] 



(U + S)^ 



-S,; I Voi P \ 5o(0, rrii, Mw) - BoiMw, rrii, M, 



w) 



M'w 



-^22(0, mi, Mw) + B22{Mw, mi, Mw) 



(7) 



The Standard Model contributions to the gauge boson 2— point functions can be found from Appendix 1 
of Ref. [27! by setting (5 = 7 = and dropping the contributions involving _fC° and . Our convention 
in this note for the sign of the 2-point functions is opposite from that of Ref. [27|] . 



4 



The definitions of tlie Passarino-Veltman B functions are given in Appendix A. The contri- 
butions from the Goldstone bosons are identical to the Standard Model case and hence are 
not included in Eq. [7f . 

Using the results in Appendix A, 

9 - V 11/ |2 / I \ ]r^^^'^ 



+Fi(M|, m,, Mz) - J]^F2{Ml m„ Mz) + Csj 



^u+s),- -y.\v,,\ [-imt^ Mf-M^^V-mi) 

i2 



2 



+Fi(M^, m,, Mw) - —j^F^iM^, rrii, Mw) + Cu\ , (8) 

where the terms Cs, Ct and Cu represent contributions which are independent of rrii. 

In order to compare with fits to data, we must subtract from Eq. [7] the Standard Model 
Higgs boson contribution evaluated at a reference Higgs mass, Mf{,ref, 



TT 



""'"'^ -1 « Ml ^Ml^.,-Ml\ AMI J \ Ml 



+Fi(M|, MH,.e/, Mz) - -^^F2{Ml, MH,ref, Mz) + Cs 



...... _ If l ^^,.e/ , M^,,,^ / Ml^^f^ f 

+Fi(M^, MH,.e/, M^) - -^F2(M^, MH,ref, Mw) + Cu\. (9) 
Finally, we compare the quantities from Eqs. [7] and [9], 





— •S'^ 


— SH,ref 






= T, 


~ TH,ref 






= 


— UH,ref , 


(10) 



with a fit to experimental data in order to obtain limits on the allowed masses and mixing 
angles. For rrii, MH,ref » Mw, Mz, we find the famihar forms ^24]. 

AS. = —Ei I K)^ I' log ^ 



12vr-" —\Ml,^, 
^'^<t> = ~T7 — I ^Oi 1^ log 



AU^ = 0. (11) 



Eq. [7] is in agreement with the results of Ref. when the Goldstone boson contributions are included. 
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For rrii ~ Mw,Mz the 0{—^,—§-) terms which are neglected in Eq. [11] are numerically 

i i 

important^. Our fitting proceedure includes the complete result and is described in Appendix 
B. 



IV. RESULTS 

In this section, we consider models with one and two scalar singlets in addition to the 
SU (2)l doublet, and extract the regions of parameter space allowed by precision electroweak 
measurements. The goal is to draw some general conclusions about the Higgs discovery 
potential in models with expanded scalar sectors. 

The dominant discovery channel for much of the Higgs mass range is 0j ZZ* — ^ 4 
leptons. The production rates of the (pi are reduced from the Standard Model rates by 
I Voi p. For irii < 200 GeV , the 0^ scalar decay width is less than or comparable to the 



detector resolution |29l. |30[ |. so we use the narrow width approximation and neglect effects 
of the finite scalar widths. For the lightest Higgs boson, the Higgs branching ratios are 
identical to the Standard Model branching ratios. For the heavier Higgs bosons, the scalar 
branching ratios depend on whether the 0j — > cpjCpk channel is accessible for some (f)j and 
(pk- Whether or not this channel is open depends on the scalar mass spectrum, along with 
the parameters of the scalar potential. We define Qjk = 1(0) if the decay (pi (pj(pk is (is 
not) allowed. The signal for (pi production with the subsequent decay to Standard Model 
particles is then suppressed from the Standard Model rate by[2i], 

IT/ 12 p5M 

' ' IK). prf^^ + s,,o,fcr(0,-.0,0,)' ^'^^ 

where F^*^ is the total width in the Standard Model for a Higgs boson of mass m^. For 
Cijk = 0, Xj = Voi- From Eq. [121 is always less than one, so the addition of scalar singlets 
reduces the significance of the usual Higgs discovery channels. 

Fig. [1] shows the minimum value of Xj, X™", for which a 5a significance in the (pi — >■ 
ZZ* 4 lepton channel can be found at the LHC with ^/s = 14 TeV and L = 30 fb~^. This 
figure is obtained by rescaling recent ATLAS studies[30]. As long as the (pi (pjCpk channel 
is closed for the heavier scalars, then this limit can be trivially applied for all (pi and Voi^. 
For a model with one singlet and one SU{2)l doublet (and hence 2 physical Higgs bosons), 
if both scalars have masses less than m^. ~ 160 GeV, then since at least one of the scalars 
must have Voi > l/\/2, at least one scalar can be discovered through the (pi — >■ ZZ* — > 4 
lepton channel. The situation changes when a second singlet is added. Now there are three 
physical scalars and it is possible for all scalars to have masses less than ~ 160 GeV and to 
have mixing angles V^i ~ l/"\/3- In this case none of the scalars will be seen (at least with 
L = 30 fb~^) in the (pi — > ZZ* 4 lepton channel. This is a generalization of the result of 
Ref. in and can be straight-forwardly applied to examples with more singlets. 

In the region 165 GeV ^ rrii ^ 180 GeV, the (pi —>■ ZZ* — > 4 lepton channel does not 
lead to a 5a discovery with 30 fb^^. In this mass region, the most useful discovery channel 
is (pi — i> W^W~ — >• e^vjji^u, which yields a > 5cr discovery for 140 GeV ^ rrii < 185 GeV 
for Xj ~ ifiO]. For rrii ~ 160 GeV, a 5a discovery is possible with Xj > 0.7. 



^ Ref. [if retains only the logarithmic contributions. 

^ This also requires that the mass differences between the scalars be greater than the detector resolution. 



6 



5a Discovery from (p -> ZZ -> 4 leptons 




200 250 300 

m, (GeV) 



FIG. 1: Minimum value of Xi for which a 5a significance in the (pi ZZ* 4 lepton channel is 
obtained with the ATLAS detector at the LHC with ^/s = U TeV and / L = 30 /6~^[30]. 



A. Fit with One Singlet 

Figs. [2] and [3] show results with one singlet scalar in addition to the SU (2) ^ scalar doublet. 
The scalar sector is described by the masses of the two scalars, tuq and mi, and one mixing 
angle which we take to be Vqi (Voo = y^l — V^^i). The fit to the experimental limits on AS", 
AT and AU is performed as described in Appendix B and the maximum allowed value of Vqi 
for various values of mo is shown in Fig. O (For simplicity, we assume Qjk = for all i,j, k.) 
For Vol ~ 0, 00 is predominantly the neutral component of the SU{2)l doublet with nearly 
Standard Model couplings and the 95% confidence level limit on the allowed value of mo is 
just the 95 % confidence level hmit of this fit in the Standard Model, Mh,sm ^ 166 GeV. 
There is no limit on mi in this case. 

For moderate mixing, the heavier scalar, (pi, can be quite heavy. For example, the lightest 
scalar could have mo ~ 140 GeV with a coupling Vqo ~ -7 while the heavier scalar could 
have a mass mi ~ 200 GeV with a coupling Vqi ~ .7. In this case, comparison with 
Fig. [1] shows that both scalars could be observed in the ZZ* —>■ 4 lepton channel with 
30 fb~^. If (pi becomes too heavy (say (pi ~ 500 GeV), then its coupling to Standard Model 
particles is restricted by the precision electroweak measurements to be less than Vqi ^0.5 
(for mo ^ 114 GeV) and so (pi cannot be found in the ZZ* 4 lepton mode with 10 fb~^. 

Fig. [3] demonstrates that scalars which have masses in the 200 GeV range can be com- 
patible with the electroweak precision measurements and have couplings large enough to be 
discovered at the LHC. In much of the parameter space of this plot, both scalars will be 
observed. 



7 



1 Singlet 

1 
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FIG. 2: Allowed region (at 95% confidence level) in a model with one additional singlet in addition 
to the usual SU{2)l doublet. The lightest (heavier) scalar is mo {mi) and the mixing matrix is 
defined in Eq. [H The region below the curves is allowed by fits to S, T and U . 



B. Fit with Two Singlets 

In this subsection, we examine how the allowed masses of the scalars are changed with 
the addition of two singlets in addition to the Standard Model doublet. The scalar sector 
now has three scalars with masses mo, mi and m2 and the mixing matrix ^ is a 3 x 3 unitary 
matrix. The phenomenology is quite different from the case with one singlet. As mentioned 
previously, with two singlets it is possible to sufficiently suppress the couplings V^i to all 
scalars such that none of them are observable with 10 fh~^ at the LHC if they all satisfy 
the Standard Model limit, rrii < 166 GeV . 

Figures H] and O show the minimum allowed value from the electroweak fit for Vqi as a 
function of Vqo for fixed masses. (We assume Qjk = for simplicity). The minimum of 

Vol results from requiring that the coupling to the heaviest scalar, V02 = y^l — Vqq — Vqi, 
not be large enough that 02 makes a significant contribution to AS", AT or AU. The solid 
red lines in Figs. [Hand [5] are Voi = 0.6 which roughly represents the limit of observability 
in the (pi ZZ* 4 leptons channel. In these examples, there is never more than one 
scalar is observable. In the regions enclosed by the dotted lines, all three scalars would elude 
detection in the (pi — > ZZ* — »• 4 leptons channel with 30 fb~^. 

The heaviest scalar can have a mass in the m2 ~ 200 — 250 GeV range and still have a 
coupling, V02, large enough to be observed in the ZZ* —>■ 4 lepton channel if mo and mi 
are less than ~ 160 GeV, although the lighter scalars will have couplings which are too 
small to be observed in this example. Thus observation of a scalar Higgs-like particle with 
mi > Mn^sm can be considered as a smoking gun for theories with multiple scalar singlets. 
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1 Singlet 




m„ (GeV) 

FIG. 3: Allowed region (at 95% confidence level) in a model with one additional singlet in addition 
to the usual SU{2)l doublet. The lightest (heavier) scalar is mo {jni) and the mixing matrix is 
defined in Eq. HI The region below and to the left of the curves labelled with values of Vqi are 
allowed by fits to S, T and U . The region to the right of the dashed line is allowed by direct search 
limits from LEP2. 



V. CONCLUSIONS 

We have considered the discovery potential for Higgs bosons in theories with multiple 
scalar singlets and demonstrated that quite simple modifications of the Standard Model 
Higgs sector can reduce the significance of the standard Higgs discovery channels. 

The addition of two scalar gauge singlets can change the Higgs sector dramatically from 
the Standard Model and also from the case with a single scalar. In this case it is possible 
to hide the Higgs boson if all three physical scalars are light, m-j ^ 160 GeV , with roughly 
equal mixing angles, Vqi ~ Alternatively, in this case, the electroweak precision 

measurements allow a Higgs boson in the 200 — 250 GeV mass region with couplings to 
Standard Model particles which are large enough to allow discovery. 
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Appendix A 



The Passarino-Veltman functions|3ll]. are defined as, 

;Bo{q^, nil, 1712) = J 



167r2 



(27r)" [k^ - ml][{k + - ml 
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2 Singlets 
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FIG. 4: Allowed region at 95% confidence level for a model with two singlets in addition to the 
SU{2)l scalar doublet. The allowed regions are above and to the right of the dashed and dot- 
dashed curves. The solid curve is Yl,i I^OiP = 1- The curved dotted line is V02 = 0.6, while the 
straight dotted lines are Vbo = 0.6 and Vqi = 0.6. 



167r2 



(27r)" [P -mf][(A; + g)2 -m^ 



(13) 



We define, 



Bo{p^,mi,m2) = [A^2]|^ - -^l(p^"^l,"^2)| 



12m 



2(; + l) -^^2(p',mi,m2)} (14) 



where 



Fi(p^,mi,m2) = / dx\n(l — x + — 
Jo \ r 



X p^x^l — x) 



ml 



(1 — x) 



^2(^^,7711,7722) = / dxl (1 — x) + rx — ^^x(l — x)\ ln( X + — -^x(l — x) {lb) 

JO I mf ) \ r m2 



We need the special cases [32 
Fi(0,mi,m2) = 
F,{mlm„m2) = 



1 — r 



ln(r) 



2-— log(r) + ^lnf^ ^ 

2r ^ 2r V1 + /3 
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2 Singlets 




^00 

FIG. 5: Allowed region at 95% confidence level for a model with two singlets in addition to the 
SU{2)l scalar doublet. The allowed regions are above and to the right of the dashed and dot- 
dashed curves. The solid curve is 

^. |Vb.|2 = 1. The solid curve is X^il^OiP = 1- The curved 
dotted line is Vq2 = 0.6, while the straight dotted lines are Vbo = 0.6 and Vqi = 0.6. 



F2(0,mi,m2) 



F2{ml,mi,m2] 



1+r 



2(1 - r 



ln(r) 



1 

4r 



1 + Fiiml, nil, m2 



or 



1 2r 
3 + ¥ 



(16) 



and we define 



r = 



(3 ^ 



rrin 



No 



V ml 



ra + e) 



(17) 



Appendix B 

We use the fit to electroweak precision data given in Ref . Q , 

AS = S-SsM = -0.126 ±0.096 
AT = T-TsM = -0.111 ±0.109 

A^7 = ?7- f/cjM = ±0.164 ±0.115 (18) 
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with the associated correlation matrix, 

/l.O 0.866 -0.392 \ 
Pij = 0.866 1.0 -0.588 . 
\ -0.392 -0.588 1.0 / 

Ax^ is defined as 

Ax' = S.,(AX, - AX,){a%\AX, - AX,) , (19) 

where AX, = AS, AT, and AU are the central values of the fit in Eq. [TSl AXj = AS^, AT^, 
and AUff) from Eq. (TDl (Ti are the errors given in Eq. [TS] and af^ = criPijaj. The 95% 
confidence level limit corresponds to Ax^ = 7.815. We vary the input values of Voi and 
to find the Ax^ = 7.815 contours shown in Figs. 2-5. This fit gives a 95% confidence level 
limit on the Standard Model Higgs boson of Mh^sm < 166 GeV . 
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